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1 ^ . Abstract 

. We evidence the existence of plasma oscillations of electrons-positron pairs created 

by the vacuum polarization in an uniform electric field with E < E c . Our general 
treatment, encompassing also the traditional, well studied case of E > E c , shows 
the existence in both cases of a maximum Lorentz factor acquired by electrons 
and positrons and allows determination of the a maximal length of oscillation. We 
quantitatively estimate how plasma oscillations reduce the rate of pair creation and 
increase the time scale of the pair production. These results are particularly relevant 
in view of the experimental progress in approaching the field strengths E < E c . 
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The plasma oscillation phenomenon of the electron-positron pairs created by 
vacuum polarization for E > E c = m 2 c 3 /(eh) (H is the Planck's constant, c 
is the speed of light, m and e are the electron mass and the absolute value 
of electron charge respectively) represents one of the most popular topics in 
relativistic field theory today. In particular the Vlasov-Boltzmann equation 
has been used within QCD in [1],[2], and with the semiclassical field equations 
in [3] within QED. This approach was shown in [3] to be in excellent agreement 
with quantum field theory calculations [4] . Applications of these studies range 
from heavy ion collisions [5]- [7] to lasers [8]. 

We have introduced collisional terms in the Vlasov-Boltzmann equation for 
such a system in [9]. Our results have been considered of interest in the studies 
of pair production in free electron lasers [10], [11], in optical lasers [12], of 
millicharged fermions in extensions of the standard model of particle physics 
[13], electromagnetic wave propagation in a plasma [14], as well in astrophysics 
[15]. 
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In this Letter we explore the case of undercritical electric field which has not 
yet been studied in the literature. It is usually expected that for E < E c back 
reaction of the created electrons and positrons on the external electric field can 
be neglected and electrons and positrons would move as test particles along 
electric lines of force. Here we show that this is not the case in a uniform 
unbounded field. This work is urgent since the first observation of oscillations 
effects should be first detectable in experiments for the regime E < E c , in 
view of the rapid developments in experimental techniques, see e.g. [16], [17]. 

We introduce an approach based on continuity, energy-momentum conserva- 
tion and Maxwell equations in order to account for the back reaction of the 
created pairs. By this treatment we can analyse the new case of undecritical 
field, E < E c , and recover the old results for overcritical field, E > E c . In 
particular, we are focusing on the range 0.15E C < E < 10E C . 

It is generally assumed that electrons and positrons are created at rest in pairs, 
due to vacuum polarization in uniform electric field with strength E [18]- [23], 
with the average rate per unit volume and per unit time 
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This formula is derived for uniform constant in time electric field. However, 
it still can be used for slowly time-varying electric field provided the inverse 
adiabaticity parameter [22]- [25] is much larger than one, 
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where u is the frequency of oscillations, T = m/u is dimensionless period of 
oscillations. Equation (2) implies that time variation of the electric field is 
much slower than the rate of pair production. In two specific cases considered 
in this paper, E = 10E C and E = 0.15E C we find for the first oscillation 
i] = 334 and rj = 3.1 x 10 6 respectively. This demonstrates applicability of the 
formula (1) in our case. 

From the continuity, energy-momentum conservation and Maxwell equations 



1 We use in the following the system of units where h = c = 1, e = yfa. ~ a/1/137, 
a being the fine structure constant. 
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written for electrons, positrons and electromagnetic field we have 

sf = (4) 

where n is the comoving number density of electrons, T^ u is energy-momentum 
tensor of electrons and positrons 

= run (ufaUfa + UfaUfa) , (6) 

F >iV is electromagnetic field tensor, J M is the total four-current density, W 1 is 
four velocity respectively of positrons and electrons 

Ufa = U» = 7 (1, v, 0, 0) , Ufa = 7 (1, -v, 0, 0) , (7) 

v is the average velocity of electrons, 7 = (1 — v 2 ) is relativistic Lorentz 
factor. Electrons and positrons move along the electric field lines in opposite 
directions. 

We choose a coordinate frame where pairs are created at rest. Electric field in 
this frame is directed alone; and introduce coordinate number density 

n = nj. In spatially homogeneous case from (3) we have 

n = S. (8) 

With our definitions (6) from (4) and equation of motion for positrons and 
electrons 

m-^ = TeF^ (9) 

we find 

— = -en (Ufa - Ufa) F> + mS (Ufa + Ufa) = -F^\ (10) 

where the total current density is the sum of conducting J^ ond and polarization 
Jp 0l currents [6] densities 

J* 1 = Jcond + Jpoi-> (H) 
JL* = ^ (Ufa - Ufa) , (12) 

^ = ^7(0,1,0,0). (13) 



Energy-momentum tensor in (4) and electromagnetic field tensor in (5) change 
for two reasons: 1) electrons and positrons acceleration in the electric field, 
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given by the term J^, nd , 2) particle creation, described by the term Jp d . Equa- 
tion (3) is satisfied separately for electrons and positrons. 

Defining energy density of positrons 

p = It 00 = mn 7 , (14) 

we find from (4) 

p = envE + rwyS. (15) 

Due to homogeneity of the electric field and plasma, electrons and positrons 
have the same energy and absolute value of the momentum density p, but their 
momenta have opposite directions. Our definitions also imply for velocity and 
momentum densities of electrons and positrons 

v = (16) 
P 

and 

p 2 =p 2 + m 2 n 2 , (17) 

which is just relativistic relation between the energy, momentum and mass 
densities of particles. 

Gathering together the above equations we then have the following equations 

h = S, (18) 



p = E(env + ^y (19) 

p = enE + rrnryS, (20) 

E = -8tt [env + ^) . (21) 

From (19) and (21) we obtain the energy conservation equation 



\f- + 2p = 0, (22) 

where Eq is the constant of integration, so the particle energy density vanishes 
for initial value of the electric field, Eq. 

These equations give also the maximum number of the pair density asymp- 
totically attainable consistently with the above rate equation and energy con- 
servation 

- = w (23) 

For simplicity we introduce dimensionless variables n = m 3 h, p = m 4 p, p = 
m 4 p, E = E C E, and t = m~ l t. With these variables our system of equations 
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(18)-(21) takes the form 

dh _ ~ 
4 

d t = nEv + 75, (24) 
at 

-t = nE + wS, 
dt 

dE ( >yS\ 

= — 81m nv H — — , 
dt V E J 

where S = -^E 2 exp [— -|) , v — | and 7 = (1 — v 2 )" 1 ^ 2 , a = e 2 /(hc) as 
before. 

We solve numerically the system of equations (24) with the initial conditions 
n(0) = p(0) = v(0) = 0, and the electric field E(0) = E . 

In fig. 1 we provide diagrams for electric field strength, number density, ve- 
locity and Lorentz gamma factor of electrons as functions of time, for initial 
values of the electric field E = 10E C (left column) and E = 0.15i? c (right 
column). Slowly decaying plasma oscillations develop in both cases. We es- 
timated the half-life of oscillations to be 10 3 t c for Eq = 10E C and 10 5 t c for 
E = 0.8E C respectively. The period of the fist oscillation is 50t c and 3 x 10 7 t c , 
the Lorentz factor of electrons and positrons in the first oscillation equals 75 
and 3 x 10 5 respectively for E = 10E C and E = 0.15E C . Therefore, in contrast 
to the case E > E c , for E < E c plasma oscillations develop on a much longer 
timescale, electrons and positrons reach extremely relativistic velocities. 

In fig. 2 the characteristic length of oscillations is shown together with the 
distance between the pairs at the moment of their creation. For constant elec- 
tric field the formation length for the electron-positron pairs, or the quantum 
tunnening length, is not simply m/(eE), as expected from a semi-classical 
approximation, but [26], [27] 



r>' = ^m li2 - (25, 



Thus, given initial electric field strength we define two characteristic distances: 
D* , the distance between created pairs, above which pair creation is possible, 
and the length of oscillations, D = or, above which plasma oscillations occur 
in a uniform electric field. The length of oscillations is the maximal distance 
between two turning points in the motion of electrons and positrons (see fig. 
2). From fig. 2 it is clear that D ^> D* . In the oscillation phenomena the larger 
electric field is, the larger becomes the density of pairs and therefore the back 
reaction, or the screening effect, is stronger. Thus the period of oscillations 
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Fig. 1. Electric field strength, number density of electrons, their velocity and Lorentz 
gamma factor depending on time with Eq = 10E C (left column) and Eq = 0.15-E C 
(right column). Electric field, number density and velocity of positron are measured 

respectively in terms of the critical field E c , Compton volume ij? = (^) , and the 
speed of light c. We define the length of oscillation as D = ct, where r is the time 
needed for the first half-oscillation, shown above. 

becomes shorter. Note that the frequency of oscillation is not equal to the 
plasma frequency, so it cannot be used as the measure of the latter. Notice 
that for E <C E c the length of oscillations becomes macroscopically large. 

At fig. 3 maximum Lorentz gamma factor in the first oscillation is presented 
depending on initial value of the electric field. Since in the successive oscil- 
lations the maximal value of the Lorentz factor is monotonically decreasing 
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Fig. 2. Maximum length of oscillations (black curves) together with the distance 
between electron and positron in a pair (red curve) computed from (25), depending 
on initial value of electric field strength. The solid black curve is obtained from solu- 
tions of exact equations (24), while the dotted black curve corresponds to solutions 
of approximate equation (27). 

(see fig. 1) we conclude that for every initial value of the electric field there 
exists a maximum Lorentz factor attainable by the electrons and positrons 
in the plasma. It is interesting to stress the dependence of the Lorentz fac- 
tor on initial electric field strength. The kinetic energy contribution becomes 
overwhelming in the E < E c case. On the contrary, in the case E > E c the 
electromagnetic energy of the field goes mainly into the rest mass energy of 
the pairs. 



This diagram clearly shows that never in this process the test particle approx- 
imation for the electrons and positrons motion in uniform electric field can be 
applied. Without considering back-reaction on the initial field, electrons and 
positrons moving in a uniform electric field would experience constant accel- 
eration reaching v ~ c for E = E c on the timescale t c and keep that speed 
thereafter. Therefore, the back reaction effects in a uniform field are essential 
both in the case of E > E c and E < E c . 

We compare the average rate of pair creation for two cases: when the elec- 
tric field value is constant in time (an external energy source keeps the field 
unchanged) and when it is self- regulated by equations (24). The result is rep- 
resented in fig. 4. It is clear from fig. 4 that when the back reaction effects are 
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Fig. 3. Maximum Lorentz gamma factor 7 reached at the first oscillation depending 
on initial value of the electric field strength. 

taken into account, the effective rate of the pair production is smaller than 
the corresponding rate (1) in a uniform field E . At the same time, discharge 
of the field takes much longer time. To quantify this effect we compute the 
efficiency of the pair production defined as e = n(ts)/no where ts is the time 
when pair creation with the constant rate S(E ) would stop, and no is defined 
above, see (23). For E = E c we find e = 14%, while for E = 0.3E C we have 



It is clear from the structure of the above equations that for E < E c the 
number of pairs is small, electrons and positrons are accelerated in electric 
field and the conducting current dominates. Assuming electric field to be weak 
we neglect polarization current in energy conservation (19) and in Maxwell 
equation (21). This means energy density change due to acceleration is much 
larger than the one due to pair creation, 



In this case oscillations equations (18)-(21) simplify. From (19) and (20) we 
have p = vp, and using (16) obtain v = ±1. This is the limit when rest mass 
energy is much smaller than the kinetic energy, 7 ^> 1. 

One may therefore use only the first and the last equations from the above set. 
Taking time derivative of the Maxwell equation we arrive to a single second 



e = 1%. 



Eenv ^> rwyS. 



(26) 
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Fig. 4. The average rate of pair production n/t is shown as function of time (thick 
curve), comparing to its initial value S(Eq) (thin line) for Eq = E c . The dashed line 
marks the time when the energy of electric field would have exhaused if the rate 
kept constant. 



order differential equation 
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Equation (27) is integrated numerically to find the length of oscillations shown 
in fig. 2 for E < E c . Notice that condition (26) means ultrarelativistic approx- 
imation for electrons and positrons, so that although according to (18) there 
is creation of pairs with rest mass 2m for each pair, the corresponding increase 
of plasma energy is neglected, as can be seen from (26). 

Now we turn to qualitative properties of the system (18)-(21). These nonlinear 
ordinary differential equations describe certain dynamical system which can 
be studied by using methods of qualitative analysis of dynamical systems. 
The presence of the two integrals (17) and (22) allows reduction of the system 
to two dimensions. It is useful to work with the variables v and E. In these 
variables we have 



dv 

dt 
dE 

~dl 
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(28) 
(29) 



9 



Introducing the new time variable r 



dr 




1/2 



(30) 



we arrive at 



dv 
dr 
dE 




(31) 



(32) 



Clearly the phase space is bounded by the two curves v — ±1. Moreover, 
physical requirement p > leads to existence of two other bounds E = ±Eq. 



focus at E — and v — 0. 

The phase portrait of the dynamical system (31), (32) is represented at fig. 5. 
Thus, every phase trajectory tends asymptotically to the only singular point 
at E — and v — 0. This means oscillations stop only when electric field 
vanishes. At that point clearly 



is valid, i.e. all the energy in the system transforms just to the rest mass of 
the pairs. 

In order to illustrate details of the phase trajectories shown at fig. 5 we plot 
only 1.5 cycles at fig. 6. One can see that the deviation from closed curves 
shown by dashed curves is maximal when the field peaks, namely when the 
pair production rate is maximal. 

The above treatment has been done by considering uniquely back reaction 
of the electron-positron pairs on the external uniform electric field. The only 
source of damping of the oscillations is pair production, i.e. creation of mass. 
As our analysis shows the damping in this case is exponentially weak. How- 
ever, since electrons and positrons are strongly accelerated in electric field the 
bremsstrahlung radiation may give significant contribution to the damping of 
oscillations and further reduce the pair creation rate. Therefore, the effective 
rate shown in fig. 4 will represent an upper limit. In order to estimate the 
effect of bremsstrahlung we recall the classical formula for the radiation loss 
in electric field 



This system has only one singular point in the physical region, of the type 
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Fig. 5. Phase portrait of the two-dimensional dynamical system (31), (32). Tildes 
are ommitted. Notice that phase trajectories are not closed curves and with each 
cycle they approach the point with E = and v = 0. 

Thus the equations (19) and (20), generalized for bremmstrahlung, are 



p = E env + 
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'-e 4 mE 2 , 



p = enE + mvyS e 4 mE 2 v. 
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(35) 
(36) 



while equations (18) and (21) remain unchanged. Assuming that new terms 
are small, relations (17) and (22) are still approximately satisfied. 

Now damping of the oscillations is caused by two terms: 
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Fig. 6. Phase trajectory for 1.5 cycles (thick curve) compared with solutions where 
the Schwinger pair production is switched off (dashed curves). 

We integrate the modified system of equations, taking into account radiation 
loss, starting with E = 10E C . We present the results in fig. 7 where the sum 
of the energy of electric field and electrons-positrons pairs normalized to the 
initial energy is shown as a function of time. The energy loss reaches 20 percent 
for 400 Compton times. Thus the effect of bremmstrahlung is as important 
as the effect of collisions considered in [9] for E > E c , leading to comparable 
energy loss for pairs on the same timescale. For E < E c we expect that the 
damping due to bremmstrahlung dominates, but the correct description in 
this case requires Vlasov-Boltzmann treatment [28]. 

The damping of the plasma oscillations due to electron-positron annihilation 
into photons has been addressed in [9]. There it was found that the system 
evolves towards an electron-positron-photon plasma reaching energy equipar- 
tition. Such a system undergoes self-acceleration process following the work 
of [29]. 
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Fig. 7. Losses of the energy due to classical bremsstrahlung radiation. The energy 
density of the system of electrons, positrons and the electric field normalized to the 
initial energy density is shown without (solid line) and with (dashed line) the effect 
of bremsstrahlung. 

We can therefore reach the following conclusions: 

• It is usually assumed that for E < E c electron-positron pairs, created by the 
vacuum polarization process, move as charged particles in external uniform 
electric field reaching arbitrary large Lorentz factors. We demonstrate the 
existence of plasma oscillations of the elecron-positron pairs also for E < E c . 
The corresponding results for E > E c are well known in the literature. For 
both cases we determine the maximum Lorentz factors 7 max reached by elec- 
trons and positrons. The length of oscillations is 10 H/(mc) for E = 10E C , 
and 10 7 h/{mc) for E = 0.15E C . We also study the asymptotic behaviour 
in time, t — > oo, of the plasma oscillations by the phase portrait technique. 

• For E > E c the vacuum polarization process transforms the electromag- 
netic energy of the field mainly in the rest mass of pairs, with moderate 
contribution to their kinetic energy: for Eq = 10E C we find 7 max = 76. For 
E < E c the kinetic energy contribution is maximized with respect to the 
rest mass of pairs: 7 max = 8 x 10 5 for E = 0.15E C . 

• In the case of oscillations the effective rate of pair production is smaller 
than the rate in uniform electric field a constant in time, and consequently, 
the discharge process lasts longer. The half-life of oscillations is 10 3 t c for 
Eq = 10E C and I0 5 t c for E = 0.8E C . We computed the efficiency of pair 
production with respect to the one in a uniform constant field. For E = 
0.3E C the efficiency is reduced to one percent, decreasing further for smaller 
initial electric field. 

All these considerations apply to a uniform electric field unbounded in space. 
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The presence of a boundary or a gradient in electric field would require the 
use of parial differential equations, in contrast to the ordinary differential 
equations used here. This topic needs further study. We also estimated the 
effect of bremsstrahlung for E > E c , and found that it represents comparable 
contribution to the damping of the plasma oscillations caused by collisions [9]. 
It is therefore clear, that the effects of oscillations introduces a new and firm 
upper limit to the rate of pair production which would be further reduced if 
one takes into account bremsstrahlung, collisions and boundary effects. 

We thank the anonymous referee for helpful suggestions on the presentation 
of our results. 
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